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Abstract 

We apply the recently introduced method of hermitization to study in the 
large N limit non-hermitean random matrices that are drawn from a large 
class of circularly symmetric non-Gaussian probability distributions, thus ex- 
tending the recent Gaussian non-hermitean literature. We develop the general 
formalism for calculating the Green's function and averaged density of eigenval- 
ues, which may be thought of as the non-hermitean analog of the method due 
to Brezin, Itzykson, Parisi and Zuber for analyzing hermitean non-Gaussian 
random matrices. We obtain an explicit algebraic equation for the integrated 
density of eigenvalues. A somewhat surprising result of that equation is that 
the shape of the eigenvalue distribution in the complex plane is either a disk 
or an annulus. As a concrete example, we analyze the quartic ensemble and 
study the phase transition from a disk shaped eigenvalue distribution to an 
annular distribution. Finally, we apply the method of hermitization to develop 
the addition formalism for free non-hermitean random variables. We use this 
formalism to state and prove a non-abelian non-hermitean version of the central 
limit theorem. 
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1 Introduction 



There has been considerable interest in random non-hermitean matrices recently. Pos- 
sible applications range over several areas of physics[0, One difficulty is that the 
eigenvalues of non-hermitean matrices invade the complex plane, and consequently, 
various methods developed over the years to deal with random hermitean matrices 
are no longer applicable, as these methods typically all involve exploiting the pow- 
erful constraints of analytic function theory. (See in particular the paper by Brezin, 
Itzykson, Parisi, and Zuber|^.) In a recent paper [Q, we proposed a "method of 
hermitization" , whereby a problem involving random non-hermitean matrices can be 
reduced to a problem involving random hermitean matrices, to which various stan- 
dard methods (such as the diagrammatic method^, or the "renormalization group" 
method]^, 0, |, ^) can be applied. 

To our knowledge, the literature on random non-hermitean matricesj^, Q has fo- 
cussed exclusively on Gaussian randomness. For instance, it has been known for over 
thirty years, from the work of Ginibre|T^, that for the Gaussian probability distri- 
bution P(0) = {1/Z)exp {—NiKp'^cf)) (here, as in the rest of this paper, denotes an 
N X N random matrix with the limit — » oo understood), the density of eigenvalues 
of is uniformly distributed over a disk of radius 1 in the complex plane. In this 
paper, we point out that using the method of hermitization we can determine the 
density of eigenvalues of probability distribution of the form 

Zj 

where V is an arbitrary polynomial of its argument. Indeed, by a simple trick, we 
show that we can obtain the desired density of eigenvalues with a minimal amount of 
work, by judiciously exploiting the existing literature on random hermitean matrices. 

Due to the symmetry of -P(0) under the transformation — e*"0, the density 
of eigenvalues is obviously rotational invariant. We find that the class of probability 
distributions of the form (|1 . 1|) exhibits a universal behavior in the sense that whatever 
the polynomial V is, the shape of the eigenvalue distribution in the complex plane is 
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always either a disk or an annulus. 

In a certain sense, our work may be thought of as the analog of the work of 
Brezin et al. for random hermitean matrices [Q; they showed how the density of 
eigenvalues of hermitean matrices taken from the probability distribution P[}p) = 
{l/Z)exp[—NtTV{!f)] with V an arbitrary polynomial can be determined, and not 
just for the Gaussian case studied by Wigner and others in which V = {l/2)tnp'^. 
An important simplifying feature of the analysis in |^ is that P{f) depends only on 
the eigenvalues of if, and not on the unitary matrix that diagonalizes it. In contrast, 
the probability distribution ( |1 . 1|) for non-hermitean matrices depends explicitly on 
the GL{N) matrix S used to diagonalize (p = S^-^AS, and 5* does not decouple. 
Remarkably however, for the Gaussian P{4>), Ginibre [|10| managed to integrate over 
S explicitly and derived an explicit expression for the probability distribution of the 
eigenvalues of 0. Unfortunately, it is not clear how to integrate over S and derive the 
expression for the eigenvalue probability distribution for non-Gaussian distributions 



of the form ( |1 . 1| ) . In this paper we circumvent this difficulty by using the method of 
hermitization. 

As an explicit application of our method, we study the case = m?(f)^(l) + 

^{(f^(f)Y in detail. For m? positive, we expect we would get a disk-like distribution 
generalizing Ginibre's work. As we make im? = —fi"^ more and more negative, we 
expect a phase transition at some critical value /i^, after which we might imagine the 
disk fragmenting into an annulus. Indeed, we find the critical value 

We also calculate the density of eigenvalues inside the annulus in detail. 

The problem of adding random hermitean matrices has been much discussed in the 
recent literature |T2|, |^, |T^, |T^, |T6[. In particular, a Feynman diagrammatic proof 
of the formalism, written in terms of the so-called Blue's function, was given in . 
In a series of very interesting papers]^, Zahed and his collaborators have extended 
the addition formalism to random non-hermitean matrices. It turns out that the 
formalism we used above to determine the density of eigenvalues of random non- 
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hermitean matrices can be naturally used to study the addition formalism. Thus, we 
obtain, in a way which we naturally think is quite transparent, the addition formalism 
for random non-hermitean matrices, thus reproducing the result of Zahed et al. We 
use this addition formalism to state and prove a non-abelian non-hermitean version 
of the central limit theorem. 
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2 The Two Dimensional Gas of Eigenvalues 



As is well known, the Dyson gas approach represents one of the most powerful methods 
in analyzing hermitean random matrices. Diagonalizing a hermitean matrix cj) = 
S'^AS one writes 



jd(t>f{ct>) = j dfi{s) fudh\m^^-^jy] fi^k,s). (2.1) 

For non-hermitean random matrices, one can diagonalize = S^^AS and write 

/(A,, A*, 5). (2.2) 



/ #d(/.t/(0,0t) = f dfi{S)Y[dX,dXl (n |A. - A,n 



The fourth power in (Aj — Xj) here, in contrast to the second power in (|2.1|) , is easily 
understood by dimensional analysis: a non-hermitean complex matrix has twice as 
many real variables as a hermitean matrix. The important difference between the 
hermitean and the non-hermitean cases is already apparent in the Gaussian case 
considered by Ginibre. The probability distribution 

P(0t0) = le-^t"^ *^<^ = le-^t-- [5StA*(5St)-A] (2.3) 
Z Z 

depends on S. In the case of hermitean matrices, S would be unitary and hence 
disappear from (|2.3|) entirely. Here, however, one must integrate over S. After some 
work, Ginibre IllOi showed that 



/ dSe-"" tr[55tA*(5St)-A] ^ ^^^^^^^^ -Q |;^^ _ A^r^e'^S. lA.P (2.4) 

and was thus able to proceed with the gas approach. Once again, the problem reduces 
to the statistical mechanics of a gas, this time in two dimensions, with a logarithmic 
repulsion between the gas molecules and confined by a harmonic potential. Following 
Brezin et al. 0] in evaluating the integral by steepest descent, we obtain 

2 . 



x= d y p{y) ^ 

\x-y\ 
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Making an analogy with two dimensional electrostatics mentioned, the solution of 
this equation is immediate: p is equal to l/vr inside a disk of radius 1 centered at the 
origin, and vanishes outside the disk. 

We see however, that had we been faced with a more complicated probability 
distribution, of the form P(0) = ^e"^ ^('^'<^) = ^e"^ ^[55tA*(sst)-iA]^ 
know how to carry out the integration over 5*. In this paper we circumvent this 
difficulty and present an alternative way of calculating the density of eigenvalues of 
0, based on the method of hermit izat ion. It is an interesting challenge to obtain our 
results presented below using the gas method, and we conjecture to that end, that in 
the non-Gaussian case there should be a closed form expression for the integral over 
S that is analogous to ( |2.4| ). 
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3 The Method of Hermitization 



Here we briefly review the discussion in our recent paper The averaged density of 
eigenvalues 

= E ^(^ - Re A,) 5{y - Im A,)) (3.1) 

i 

of the non-hermitean matrix 0, may be determined from the following two alternative 
quantities. The flrst quantity is 

F{z,z*) = (Itr log {z - <P){z* - 0t)) , (3.2) 

in terms of whichQ 

p{x,y) = -dd* F{z,z*). (3.3) 

TT 

The other quantity is the Green's function associated with (j), namely 
in terms of which 

p{x,y) = -d* Giz,z*). (3.5) 

TT 

The probability distributions (|LT|) studied in this paper are invariant under 
e*°0, rendering 

p{x,y) = p{r)/ 27V (3.6) 

circularly invariant. Rotational invariance thus leads to a simpler form of the deflning 
formula (|3.4|) for G{z, z*) which reads 

r 

7(r) = zG{z, z*) = J r'dr' p(r') , (3.7) 



whence 

p(r) = -§^. (3.8) 
r ar 



^We use the following notational conventions: for z — x + iy we define d = — ^ — i-^ 

so that dz = 1. Similarly, we define d* = = i + i-^^ , so that d*z* — 1 and also d*{l/z) = 
7r(5(2)(z). Finally, we denote \z\ = r. 
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Clearly, the quantity 7(r) is a positive monotonically increasing function, which sat- 
isfies the obvious "sum-rules" 



7(0) = and 7(00) = 1 . (3.9) 

In particular, observe that the first condition in ( |3.g| ) insures that no 6{x)6{y) spike 
arises in p{x,y) when calculating it from (|3.5| ) with G{z,z*) given by (|3.71 ), as it 
should be. 

We have mentioned in the introduction that a straightforward diagrammatic 
method is not allowed for non-hermitean matrices. The method of hermitization^ 
enables us to arrive at a diagrammatic method indirectly by reducing the problem 
of determining the eigenvalue density of random non-hermitean matrices to the more 
elementary problem of determining the eigenvalue density of random hermitean ma- 
trices (for which the diagrammatic method may be applied.) To make the paper 
self-contained, we summarize the method of hermitization here in the form of an 
algorithm. 

Given a set of random non-hermitean N x N matrices 0, construct the auxiliary 
set of random hermitean 2N x 2N matrices 

/ <j)-z\ 
H=\ . (3.10) 

-z* J 

Then one applies one's favorite method of hermitean random matrix theory to calcu- 
late the propagator associated with H, namely, 

^."(^;^,^*) = ((-^)'), (3.11) 

where is a complex variable and the indices /i and u run over all possible 2N values. 
In particular, in the diagrammatic method, one simply expands Qj^ in powers of l/rj, 
with interaction vertices H. This is a well defined procedure for large 77, and it 
converges to a function which is analytic in the complex 77 plane, except for the cut 
(or cuts) along the real axis which contain the eigenvalues of H. After summing this 



7 



series (and thus determining G^{ri; z, z*) in closed form), we are allowed to set r/ ^ 
in 

As mentioned, we can now calculate p{x, y) by two different methods, each with 
its advantages. The first is to simply observe that 



1 

H 



( 



from which it follows that 
G{z,z*) 



tr 












6^(0; 2, 



(3.12) 



enabling us to use ( p.5| ) to express p(x, y). This observation is the basis for many of 
the calculations in 0. 

An alternative is to take the trace of Q^: 



Giv;z,z*) = 7^(tr(^^j 



1 

2N 



:(tr. 



1 



(3.13) 



ri-H' AT^ w r]^ - (^z* - (t)^){z - 

from which one can construct an integral representation Q of (p.2|) , thus enabling us to 
use ( p.3|) to determine p{x, y). The determination of Q{rj] z, z*) involves only the well- 
known and more elementary problem of determining the density of eigenvalues of the 
hermitean random matrices which is the essence of our method of hermitization. 
Observe from (|3.13|) that 



^(0 - zO; z,z*) = ^(tr,^, 5 [^^{z* - 0t)(;, - 0))) (3.14) 
and thus counts the (average) number of zero-eigenvalues of the positive semi-definite 



hermitean matrix y {z* — (j)^) (2; — 0) fj] . Thus, if (on the average) has no eigenvalues 
equal to 2;, or in other words, if the density of eigenvalues of vanishes at z, then 



^(0 -iO;z,z*) = 



(3.15) 



^Speaking colloquially, we may say that the crucial maneuver here is that while we cannot expand 
in powers of z, we can arrange for z to "hitch a ride" with rj and at the end of the ride, throw 77 
away. 
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independently of the large limit. In particular, the boundaries of the eigenvalue 
distribution of in the large limit, are the curve (or curves) in the complex plane 
which separate regions where Q{0]z,z*) = from regions where Q{0]z,z*) ^ 0. As 
a matter of fact, one can infer the location of these boundaries even without an 
explicit knowledge of ^(0; z*) and G{z^ z*\ by investigating the "gap equation" for 
^(0; z, z*\ namely, the trace of both sides of ( |3.11| ), and then setting r/ = 0. The "gap 
equation" is a polynomial in Q with real coefficients (which depend on z explicitly 
as well as through G(z, z*), and the parameters of the potential l^.) Due to the 
chiral nature of if, this polynomial contains a trivial factor of ^(0;^;, 2*) which we 
immediately factor out. Setting ^ = in the remaining factor we obtain an equation 
for the boundary. This is discussed in detail in 0]. 

To summarize, the method of hermitization allows us to reduce the problem of 
dealing with random non-hermitean matrices to the well-studied problem of dealing 
with random hermitean matrices. Given the non-hermitean matrix 0, we study the 
hermitean matrix H instead. By whatever method one prefers, once one has deter- 
mined the quark propagator Q'^ (or its trace, the Green's function ^(77; 2;*)), one 
can in principle obtain p{x,y). Whether that can be done in practice is of course 
another story. 

It is worth emphasizing that the formalism described here has nothing to do 
with large A^ as such. The formalism is also totally independent of the form of the 
probability distribution P{(f),(f)^). 
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4 Non-Gaussian Ensembles 



Our task in this paper is to solve the more problem of determining the Green's function 
and the averaged eigenvalue density of a non-hermitean random matrix 0. Here we 
show that by applying a simple trick, we can obtain the desired density of eigenvalues 
with a minimal amount of work, by judiciously exploiting the existing literature on 
random hermitean matrices. 

In this section we will follow a common practice and borrow some terminology 
from quantum chromodynamics: we may consider 0, 0"!" as "gluons" (in zero space- 
time dimensions), which interact with a 2N dimensional multiplet of "quarks" if)^, 
with a complex mass matrix (the "inverse propagator" ) 

( Tj Z \ 



Go' 



(4.1) 



\ Z* T] ) 

(expressed in terms of its x blocks.) The matrix Green's function 



{{1 / {r] ~ H))'^) , in (|3.11|) , which represents the quark propagator, is given by the 



Dyson-Schwinger equation in terms of the one-quark irreducible self energy S(^: 

We emphasize that this equation always holds, and practically amounts to the 
definition of the one-particle irreducible self energy. For probability distributions of 
the form P(0) = (i/^)e-^t=-^(<^^<^) in it is easy to see thatQ 

Sj: = S(r7;r)5^: (4.3) 

(where we recall that r = \z\.) Note that S depends on r and not separately on z 
and z* because it may be expanded as an infinite sum of traces of Qq, each term of 



^Split the 2N dimensional quark muhiplet into two N dimensional "flavors", ip = {u,d). Then 

note that the "quark-quark-gluon" interaction is ^/'^ ^ ^| ) ^ ~ u^(l)d + d^(f)^u, and use the fact 

that the fermion lines are directed, e.g. a (j) vertex always absorbs a d and emits a u. Then invoke 
"index democracy" . 
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which is a function only of rj and r. Thus, 



s^^ i 



— {rj — S)^ 



(4.4) 



V 



leading to 



g = —yg^ = — — . (4.5) 

The important point is that in the large N limit, the one-quark irreducible self 
energy S can be written in terms of the cumulants of -P(0) = (1/Z)e~^*''^'''^^'^-* (Eq. 
( |1.1|) ), namely, the connected correlators involving k 0's and 0^'s. Diagrammatically 
speaking, the may be thought of as a "blob" out of which emanate k 0's and 0^'s, 
which cannot be separated into two smaller blobs, with ki 0's and 0^'s in one blob 
and k2 0's and 0^'s in the other (with ki + k2 = k of course.) We then obtain 

00 

J:{r^;r) = j:r2k[g{v;r)r'\ (4.6) 

k=l 

in a fashion similar to the calculation infl^. Recall also that for the Gaussian en- 
semble only the variance r2 7^ 0, and thus simply E = r2^. 

If we knew E in terms of Q, then we can solve for Q and hence the density of 
eigenvalues. We are thus faced with the problem of determining the r2A,.'s which 
appears to be a rather difficult task. 

It is important to note that the r2fc in ( [4.6| ) depend only on the probability dis- 
tribution P(0) shown above and not on the particular quantity we average over. 
Consider the problem of determining the eigenvalue density of the hermitean matrix 
We would study the Green's function 

F{w) = (-tr,,^ ——r-) = / , (4.7) 







w — a 



where p(/i) is the averaged eigenvalue density of We see from ( p.l3|) that this 
Green's function is related to ^(r/; r = 0) simply through 

g(r];r = 0) = r]F{ri'^) . (4.8) 
11 



Also, for 2; = we have from ( f4.5| ) 

j:{v;0) = v-G-\v;0)- (4.9) 

The crucial observation is that F{w) is already known in the literature on chiral 



and rectangular block random hermitean matriceqj for the Gaussian distribution |T^, 
|l^, as well as for non- Gaussian probability distributions of the form ( |1 . 1|) with 
an arbitrary polynomial potential y(0"''0))|[l^, |2^. In fact, the authors of [0, ^ 
simply calculated the diagonal elements of the propagator Ql^{rj;0) using Dyson gas 
techniques (so the coefficients are only implicit in these papers.) 

The contents of ( |4.6|) , ([4.8| ) and ( [4.9|) may thus be summarized as follows: there 
is a unique function a(^) which behaves like 1/^ as ^ tends to infinity, is regular at 
^ = 0, and satisfies the equation ^ — a^^ = T2ka'^''~^. This function is 

a{O = ^F{e)^gi^;0). (4.10) 

We, on the other hand, are interested in the opposite case, where 1] = 0,r = \z\ 
0. Consider the matrix g^{0 — iO; z, z*). Due to " color index democracy", its x 
blocks are all proportional to the unit matrix Ijy, and we see that 

( g G*\ 

gi^{0-tO;z,z*)=\ , (4.11) 

\G g I 

where ^(0 — iO; r) is given in (|3.14|) and the complex Green's function G{z,z*) was 
defined Q in ( |3.12D . Inverting the propagator in ( [4.11| ) and using ( ^?T| ) (with r/ = 0) 
and ([4.2|) , we obtain the two equations 

and 

^Namely, matrices H with the block structure i? = ^ ) 

^From this point on, unless otherwise stated, our notation will be such that Q = Cy(0 — iO;r), 
S = l](0;r) and G = G(z,z*). Note in particular from (3.14) that Q is then pure imaginary, with 
Im g > 0. 
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Note from ( [4.12| ) that if ^ = 0, then G = 1/z, which corresponds to the region 
outside the domain of the eigenvalue distribution p|, ^. In addition, using ( [4.12|) twice, 
we obtain 

\G\' 



zG 



{\G\'-g') r\ 



(4.14) 



\G\^-g^ 

Thus, due to the fact that Q is pure imaginary (see ( p.l4| )) we conclude that 7 = zG 
is always real and non-negative. Alternatively, we can write 



,2 zG {zG - 1) 



(4.15) 



and therefore from the fact that < we conclude that < 7 < 1, in accordance 



with the monotonicity of 7(r) and the sum-rules 7(0) = and 7(c)o) = 1 in (|3.9|) . 
Let us now define the pure imaginary quantity 



(4.16) 



k=l 



where the last equality follows from ( [4.6[ ) (with rj = 0). Using ( [4.12|) and ( [l.ltj| ) we 
thus have simply 

(4-17) 



It then follows from (^4.101) and the statement preceding that equation, and from 
( CT ) that 



me 



zG 
T 



(4.18) 



We have thus eliminated the ! 

From (|4.15| ) and from ( |4.17|) we determine ^ = 7/^ to be 



7 2 

r . 



7-1 



(4.19) 



Note that < 0, consistent with the bounds < 7 < 1. Comparing ( ^.18| ) and ( [4.8| ) 
we obtain the remarkable relation 



for and r that are related by ([4.19|) (with ^ imaginary.) 
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(4.20) 



We now have the desired equation for 7 = zG\ substituting (|4.19| ) into ([4.1 8|) we 
obtain 



7 



F 



7 r 



-7 + 1 







(4.21) 



Thus, given F we can solve for 7(r) using this equation. 



From (|4.71 ) we observe that for w large we may expand 



F{w) 



(4.22) 



w 

where f{w) = fo + {fi/w) + ■ ■ ■. Substituting 10 = ^"^ = 7^^/(7 — 1) into ([4.22|) and 
letting 7 tend to 1, ( [4.2 1| ) becomes 

7(^) 



(7-1) 



7r^ 



(4.23) 



and thus 7 = 1, which corresponds to the region outside the eigenvalue distribution, 
is always a double root of ([4.21| ). 

We will now explain how the Green's function F{w) is obtained in the literature]^, 
pjj , |20| . Let \^ be a polynomial of degree p. From the definition of F{w) in ([4.7| ) and 
its analyticity property, we expect F, following the arguments of Brezin et al, to have 
the formal 

F{w) = lv'{w) - P{w)^{w-a){w-b), (4.24) 



where 



P(w) 



p-2 
k=~l 



(4.25) 



The constants < a < 6 and are then determined completely by the requirement 

that F{w) ^ ^ a.s w tends to infinity, and by the condition that F{w) has at most 

an integrable singularity as w — >■ 0. [] 

^Here we assume for simplicity that the eigenvalues of (ji^tj) condense into a single segment [a, b]. 
^Thus, if a > 0, inevitably a_i = 0. However, if a = 0, then a_i will be determined by the first 
condition. 
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5 Phase Transitions in the Quartic Ensemble 



Having determined F{w) in this way, we substitute it into (^4.21|) and find G{z,z*). 
We can thus calculate the density of eigenvalues explicitly for an arbitrary V. For 
the sake of simplicity we now focus on the quartic potential^ 

V(0t0) = 2mVV + • (5.1) 

5.1 The Disk Phase 

For not too negative, which as we shall see translates here into > — y^2^, we 
expect the density of eigenvalues to be a disk. In that case we expect that p(0) > 0, 
as in Ginibre's case. Thus, from ( p.8|) p(r) = {l/r){d'~f /dr) and from the first sum 
rule 7(0) = in (|3.9|) we conclude that 07(r) ~ near r = 0. Therefore, for r small, 
( |4.21| ) yields F{—r^) ~ — 1/t^, namely, F{w) ~ 1/v^ w ~ 0, as we could have 



anticipated from Ginibre's case. This means that a = in (|4.24| ), and therefore, 
F(w) 



+ gw — (^^^—^ + ao^ \jw{w — h) . (5.2) 



From the asymptotic behavior of F{w) we then obtain]^ 

ao= g, c = a_i = , and b = . (5.3) 

3 3g 

For later convenience we record here that 

oc - I . \pA) 

Substituting (§3) and (^ into (|]^) we arriveQ at 

7^ - 7 [1 + 2(mV2 + gr% + focV^ = (5.5) 



^In the Gaussian case, V = <?!)''' </i, we have 2^i^(^'^) = ^ — ■\/^-^ — 4, whence the roots of ( 4.2l| ) are 
7 = 0, 1 and . We note that 7 = is unphysical, 7 = (i.e., G = z*) corresponds to Ginibre's 
disk||l^, and 7 = 1 is the solution outside the disk. 

^Alternatively, we can reach the same conclusion about 7(r) for r small by looking at . Ac- 
cording to (3.15), must not vanish in some neighborhood of the origin. Neither can it diverge 



there. It thus follows from (4.15) that 7(r) ~ r near r = 0. 

is positive by definition. The expression for h in ( pT^ ) is the positive root (also for w? < 0) of 
the quadratic equation which determines b. b then determines c via c = rn^ + (gb/2). 



^^In deriving (5_^) we divided through by the trivial factor (7 — 1)'^ alluded to in ( 4.23 ) 
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with roots 



1 + 2(mV^ + gr^) ± J[l + 2{m^r^ + gr'^)]^ - Abc^r^ 



(5.6) 



Let us assume that > —^J2g. Under this condition the the discriminant of ( |5.5| ) 
is non-negative, and contains a factor (r^ — 6/4)^. The root of ( |5.5| ) which satisfies 
the condition 7(0) = (see ( |3.9D ) is 7-(r) (whereas 7+(0) = 1.) However, a subtlety 
arises, since 7-(r) is monotonically increasing (as it should be, by definition) only for 
< 6/4. The root 7+(r) increases monotonically only for > 6/4. Due to the factor 
(r^ — 6/4)^ in the discriminant of ( |5.5| ), we observe that 7l(v^/2— ) = 7^(v^/2+), 
and thus 7+ and 7_ may be glued smoothly into a single monotonically increasing 
function 7(r)di^fc as 

7_(r) , < 6/4 
7+(r) , > 6/4 . 



(5.7) 



, , 2 o 2 r 2M ^c2-(m2 + 2^r2)[l + 2(mV + ^r^)] , ^, 
7rp(x, y) = m^ + 2gr^ + [sgn (- - r^)] ^ ^ )^ . (5.8) 



The eigenvalue density is then given by ( |3.8| ) as 



'[1 + 2(m2r2 + ^r4)]2 - 46c2r2 

In particular, note that p(0, 0) = hc^/Ti > 0. 

We now determine the radius of the disk rj,, namely, the boundary of the eigenvalue 
distribution. In this case there is only one boundary, and thus the r dependent 
solution of ( ^.21| ) increases monotonically from 7(0) = to 7(rb) = 1, since all the 
eigenvalues are included within r < r^. Thus, setting 7(r) = 1 in ( |5.5|) and using 
(|5.4|) we obtain 

(m^ + 6(7)^/2 _ m^im^ + 9g) 



2m^)/2g = 2- 



27g 



(5.9) 



It is possible to show that r'^ > 6/4. We thus find that the density at the boundary is 



2ri + 



2m2 



9 



-1 



(5.10) 



which can be shown to be always positive. 
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In the hmit g —>■ (with > 0, of course) c tends to and both and 6/4 
tend to l/(2m^). Therefore, i f).7\) simphfies into 

f \ _ / 2mV2 , < l/(2m2) 



Thus, from (|3.8| ) (or by taking g ^ in ( ^.8|) ) we recover Ginibre's uniform density 
p{x,y) = 2m? /tt inside a disk of radius l/(2m^), consistent with the g ^ Q hmit of 

(13). 

In the strong coupling limit g/m^ >> 1 we observe that c ~ w 2(7/3 and 5 ~ 



8/3(7. Therefore, in this limit the disk shrinks like rl ~ -y/8/27(7, which means 
that the density should blow up like 1/r^. We find indeed that in this limit ~ 
(12/57r)V6^. We also find that in this limit p(0, 0) ~ (2/7r)(2/3)=^/2^ ~ (5/27)p6 < 
Pfe, and the density grows as r increases. 

To summarize, for m? > — a/2^, the 7(r) corresponding to the ensemble (|5.1|) in 
the disk phase is thus given by ( [5.7| ) for r < and is equal to one for r > Vh. 

When becomes too negative, that is, for < — a/2^, the discriminant of (|5.5|) 
becomes negative for some range of values of r, rendering 7±(r) complex, which is 
unacceptable. This signals that the disk phase becomes unstable, and the systems 
undergoes a phase transition, as we discuss in the next section. The boundary of the 
disk phase in the m?—g plane is thus given by = —^f2g. Consider approaching this 
boundary from within the disk phase by setting vn? = ~\f2g + 5, with b positive and 
small. Then, using ( |5.3| ), we find to first order in b that c = 5/2 and h = 2^ {2/ g)—5 / g. 
In particular, at the phase boundary itself c = 0. Consequently, at the phase boundary 
simplifies into 

^^"^'^^ ~ 1 2v/2^(v/2^r2-l) , l/^<r^<rl. ^^"^^^ 



where rl = \j2/g from (|5.9|) . Thus, as we decrease S to zero, p{x,y) becomes increas- 
ingly depleted inside the disk < b{6)/4, reaching complete depletion at 6 = 0, at 
which point the disk breaks into an annulus. Finally, for later reference, let us note 
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that at the phase boundary ( |5.2| ) reads 

F{w) = -J2g + gw - g 



^w(^w-2^^. (5.12) 
5.2 The Annular Phase 

We saw that at the phase boundary = —\flg the disk configuration of the previous 



sub-section broke into an annulus. We thus expect that for rr? < —y/2g the stable 
eigenvalue distribution would be annular. 

For convenience, let us switch notations according to = — yU^, and also write 

We denote the radii of the annulus boundaries by ri < r2. In the inner void, 
r < ri, we clearly have 7(r) = 0. Thus, ([4.21|) and ( [4. 241) imply a^i^/ab = 0. The 
annulus must therefore arise for a_i = (the other possible solution a = 0, a_i ^ 
leads to a disk configuration with 7 = only at r = 0, as we discussed in the previous 
sub-section.) Therefore, 

F{w) = rn? + gw — oq \J [w — a){w — h) . (5.13) 

From the asymptotic behavior of F{w) we then obtain 

/i^ /2 , , /2 , , 

ao = g, a = and b = \- \ - ■ (5-14) 

9 \ 9 9^9 

We see that a = (2/yU^)(/i^ — fi^) which is positive for fi"^ > fi^, as it should be, by 
definition. 

We now determine 7(r) inside the annulus. Substituting ( |5.13|) and ( |5.14| ) into 



(|4.21|) and factoring out all the factors of (7 ~ 1)) we obtain a linear equation for 



7 and find 

7('^)ann«««s = 1 - 2/iV^ 25(r^ , (5.15) 
which must correspond to the region inside the annulus. 



^There is one extra such factor as compared to the disk phase. 
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The 7(r) configuration ( |5.15| ) has to match continuously to 7 = at r = ri and 
to 7 = 1 at r = r2. From the matching at r = r2 we find that 

2 /"^ 2/i2 



9 f^t 



The matching at r = ri yields a quadratic equation for rf with roots 



(5.16) 



(5.17) 



2^7 /i^ 

both of which are smaller than the expression in ( |5.16|) . In order to determine which 
of these roots corresponds to rf we recall that 7(r) must interpolate monotonically 
between its two boundary values, namely, that 7'(r) > for ri < r < r2, which is 
nothing but the statement that p(r) > 0. Substituting ( |5.15|) in ( p.8|) we find 

p{x, y) = ^ fr^ - 41 = — i^^y - 41 ' n<r<r,. (5.18) 

Note that p(r) increases monotonically. The correct normalization of ( |5.18| ) is guaran- 
teed by 7(r2) = 1, as can be checked explicitly. Finally, positivity of ( |5.18|) determines 
rf in (15.171) as 



^1 = — • (5-19) 

Note that rf = l/fil = Xj ^/^g at /i^ = p^, which coincides with the radius (squared) 
of the depleted region in the disk configuration ( |5.11|) . Also, at /i^ = /i^ we have 



r\ = 2/ = \j2/g, which coincides with the disk radius squared rl in (|5.9| ) at the 
phase boundary. Thus, at the phase boundary p^ = /i^ (|5.18|) coincides with (|5.11|) , 
namely, p(x, y) is continuous at the transition from the disk phase to the annular 
phase. 

Consider approaching the phase boundary w? = —y/2g from within the annular 
phase by setting p^ = ^/2g + 6, with 6 positive and small. Then, since all the 
expressions in ( ^.14|) are linear in p^, we find that a = 6/g and b = 2^ {2/ g) + S/g. 
Thus, at the phase boundary itself a = 0, and b = 1/a/2^. Therefore, at the phase 
boundary ( ^.13|) reads 



^(^) = -J'2g + gw-g 



\ 



w{w-2J-], (5.20) 
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which coincides with ( [5.12[ ). Thus, F{w) (and consequently, the eigenvalue density 
of 0V) is also continuous at the transition. 

It is amusing to note that the various expressions we encountered in this sub- 
section are much simpler than their counterparts in the disk phase. 



20 



6 Boundaries and Shape Universality of the Eigen- 
value Distribution: the Single Ring Theorem 

In the previous section we analyzed the probabihty distribution ( |1 . 1| ) with a quartic 



a phase transition from the disk configuration to the annulus configuration. 

This state of affairs obviously generalizes. In this class of models, the domain 
of the eigenvalue distribution is always circularly symmetric, due to the rotational 
invariance of the probability distribution (|1 . 1| ) . Consider a potential with 
several wells or minima. For deep enough wells, we expect the eigenvalues of 0V to 
"fall into the wells" . Thus, one might suppose that the eigenvalue distribution of (j) to 
be bounded by a set of concentric circles of radii < ri < r2 < ■ ■ ■ < r„^^^^, separating 
annular regions on which p(r) > from voids (annuli in which p(r) = 0.) A priori, 
it is natural to assume that the maximal number of such circular boundaries should 
grow with the degree of V , because V may then have many deep minima. 

Remarkably, however, we prove below that the number of these boundaries is two 
at the most. To reconcile this conclusion with the a priori expectation just mentioned, 
we note that while the eigenvalues of the hermitean matrix (f)^(f) may split into several 
disjoint segments along the positive real axis, this does not necessarily constrains the 
eigenvalues of (j) itself to condense into annuli. This statement is made more precise 
in the Appendix. 

We call this result the "Single Ring Theorem." To prove it, let us assume for the 
moment that there are n boundaries. It is easy to see from (|3.7| ), the defining equation 
of 7(r), that in the annular void < r < r^+i in the eigenvalue distribution, 7 is a 
constant which is equal to the fraction of eigenvalues contained inside the disk r < r^. 
Thus, the equation to determine 7 (Eq. ( [4.21| ) which we repeat here for convenience) 



must have a series of monotonically increasing constant solutions 71 < 72 < ■ ■ ■ < 1, 



potential = + We obtained circularly symmetric eigenvalue 



distributions: a disk ( p.8| ) and an annulus (|5.18 ). Making im? negative enough drove 
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which correspond to the various voids. 

In particular, from (|3.9| ) we have 7 = 1 for r > Vn^^^, namely, G = 1/z. Thus, 
7 = 1 must be a solution of ( |6.1| ), the maximal allowed constant solution. It is 
straightforward to see that 7 = 1 is indeed a solution, because we know that 



y7 — 1 J 7r2 

as 7 tends to 1. Also, if ri > 0, namely, if there is a hole at the center of the 
eigenvalue distribution, then for r < ri we must have 7 = (independently of r), 
which is obviously a solution]^ of ( |6.1| ). (On the other hand, if ri = 0, so that the 
eigenvalue distribution includes the origin, an r independent solution 7 = is of 
course only a spurious solution which should be discarded.) 

Assume now that 7 = 70 is an r independent solution of (|6.1| ). Taking the deriva- 



tive of (|6.1| ) with respect to at 7 = 70 we obtain 

,dF 



7o 



7o 



70-1 



(6.2) 



which is the condition for the existence of an r independent solution 79. Thus, there 
are two possibilities: either 70 = 0, which we already encountered, or -F(^^) -|-^^^ = 
0. This equation immediately yields 

na = ^, (6.3) 

where the integration constant is fixed by the asymptotic behavior of -F(^^) as ^ ^ 00. 
But for a generic V{(j)^(f)), F{w) is given by (Eq. (|4.7|) ) 



A ^ ' w — (b'd) J w — a 







with (|6.3| ) being the asymptotic behavior of F{w) as w = ^ 00. We thus conclude 
that ^ 00 in ( |6.3D , namely, that 70 = 1. Thus, to summarize, the only possible 
r independent solutions of ( [4.21|) are 7 = and 7 = 1, which we already discussed. 



-"^•^From the paragraph right below ( 4.24 ) it is clear that wF{w) ^ at w = 0. 
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Since no other r independent solutions arise, there can be no more than a single void 
in the eigenvalue distribution, whatever polynomial the potential V^(0V) is- The 
shape of the eigenvalue distribution is thus either a disk or an annulus. 

Interestingly enough, we can arrive at the same conclusion by invoking other 
general aspects of the method of hermitization Q], and thus providing a nice self 
consistency check of our formalism. Recall at this point that the boundaries of the 
eigenvalue distribution are given by (|3.15|) , namely, the zeros of Q. But Q is given by 
(Eq. (|4A5|) ) 

y = — ^ — , 

from which we see that at the boundaries 7 can take on only two values : zero and 
one. Since 7 is a constant in the void, by continuity, these are the only possible values 
of 7 inside any of the voids. Therefore, there may be two circular boundaries at most. 
In other words, as far as the eigenvalue density of is concerned, an ensemble of the 
form ( |1 . 1|) may have two phases at most, as we concluded earlier. 

In addition to 7 = 1 (and possibly 7 = 0), ( |6.1| ) must have other roots which 
do depend on r. Among all these other roots, we expect to find a unique root 7(r), 
which is a positive monotonically increasing function of r, that matches continuously 
at the boundaries {k < 2) to the constant roots of ( |6.1| ). The actual zG{z,z*) 
of the ensemble (|1.1|) is therefore a continuous monotonically increasing function 
of r, as required by ( p.7| ). If the eigenvalue distribution is annular with boundaries 
< ri < r2 (the case ri = corresponding to the full disk), it vanishes for < r < ri, 
rises monotonically from zero to one on ri < r < r2, and equals to one for r > r2. 
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7 The Addition Formalism for non-Hermitean Ran- 
dom Matrices 



Using the method of hermitization|^, we develop in this section the addition for- 
mahsm for non-hermitean random matrices. Our results agree with recent results 



development of this addition formalism. At the end of this section we use the addi- 
tion formalism to formulate and prove a non-hermitean non-Gaussian version of the 
central limit theorem. 

Let 01 and 02 be two NxN random non-hermitean matrices, drawn independently 
from the probability distributions Pi(0|0i) and ^2(0202)- Furthermore, let Ha {a = 
1,2) be the 2N x 2N hermitean Hamiltonians associated with (pa, as defined by (|3.10|) , 
and similarly define the "total Hamiltonian" as 



The addition problem for 0i and 02 may be formulated as follows: Consider the 
probability distributions 



of Zahed et a/. 0. We feel that the method of hermitization clearly elucidates the 




(7.1) 



^a(0l0a) = ^ exp -iVtrK(0l0a) 



a = 1,2. 



(7.2) 



a 



Given the two propagators 




(7.3) 



(where (■ ■ ■)^"'^ indicates an average with respect to Pq), find (in the large limit) 
the propagator 




(7.4) 



associated with the probability distribution 



P(0l01,0T02) = Pi(0l0i)P2(0502) . 



(7.5) 
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The propagators G^"'l^{r]; z, z*) may be obtained in principle, by applying the for- 
malism of Section 3. Once one has succeeded in solving the addition problem and 
obtained an explicit expression for ( [7.4| ), one may set r/ = and calculate the vari- 
ous averages F{z,z*) = (^tr log {z - (j)){z* - (p^)) ,G{z, z*) = {{l/N) tr l/{z - 0)) 
and Q{0;z,z*) as prescribed in Section 2, from which one may further calculate the 
averaged density of eigenvalues p{x,y) of 0i + 02. 

The formulation of the addition problem for non-hermitean random matrices 0i, 02 
bears obvious similarities to the formulation of the addition problem for hermitean 
random matrices. In fact, if in the formulation just given one replaces the bare inverse 
propagator Qq^ defined in ( [4.1|) simply by ^g^^ = z, and the chiral matrices 

/ 0, 

by the ordinary hermitean matrices ipa (a = 1, 2), one would recover the formulation 



of the usual addition problem for hermitean random matrices |13|, |16|. These iden- 
tifications of various quantities that appear in the formulations of the two addition 
problems, enable us to derive the addition formalism for random non-hermitean ma- 
trices by following the steps of the by now well known addition formalism for ordinary 



hermitean random matrices [p!4|, which we now proceed to do. 



Let T^j^ and 



^^^^ 

= j:(-)(n;r)5i^ (7.6) 

be, respectively, the connected 2k gluon "blob" (or cumulant) and the self-energy 
associated with the probability distribution Pa{4>\4>a)- From the definition of self- 
energy 



(Eq. (U)) we find 



SH'' = te' - ^"'^"T • (7-7) 
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In addition, rewriting the expansion ( [4.6|) for S^")^ we have 

-,2k-l 



Ep(a) 

k=l 



1 

2iV 



fc=i 



(7i 



For reasons that will become clear below we bothered to display the trace operation 
in ( |7.8| ) explicitly. The two gluon types 0i and 02 in the combined probability dis- 
tribution (|7.5| ) do not interact and therfore the self-energy = T,{ri;r)6^ associated 
with ( [7.5| ) is given by 

1 -I 2k— 1 oo / 



k=l 



2N 



E(Er£M[^(r];r)] 



2k-l 



(7.9) 



k=l 



where now obviously the internal quark lines are given by 

The right hand side of ( [7.8|) and of ( [7.9|) are expansions in powers of the trace 
of three different propagators: ^^""1^ (a = 1, 2) and Q^^ These three propagators 



are obviously not the same, but we observe that if they were identical, then (|7.9|) 
would simply be the sum (over the two values of a) of (|7.8|). The fact that these 
three propagators are not identical renders the relation between the three self-energies 
more complicated. The essence of the addition formalism is to maneuver around 
this difficulty. The idea is as follows: for the sake of argument, assume that it is 
mathematically consistent to consider g^"'l^{r]; ,r) as a matrix valued function of the 
inverse propagator matrix 

f V ^ \ 

\ z* r] J 



(Eq. ( [4.1|) .) This function would then have a matrix valued functional inverse Ba 
where Ba (G (Gq^^^ = Gq^- Naively, we may deduce from ( [7.7[ ) that 



Ba{G^ 



(a)aN 



-1(a) 



(7.10) 



(below we shall do away with the adverb "naively" and make the construction of Ba 
mathematically precise.) Using this matrix valued functional inverse we would be 
able to solve the equation G^^l^i^a) = Gy (where Xa is the matricial argumet of the 
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matrix valued function Q^^l^ ) as Xa = BaiG^)- We would then evaluate each of the 
T,^""^ in ( [r.8|) at the appropriate argument Xa (instead at argument ^o"^). The right 
hand side of ( [TBI ) would then be an expansion in powers of Q{ri;r), and we would 
thus conclude that Ea^^'^K'^a) = ^{Qo^)- 

The arguments in the last paragraph follow the same logic as in the case of the 
addition formalism of random hermitean matrices. Recall that in the theory of ordi- 
nary hermitean random matrices, the Blues function B{w) and the Green's function 
Q{z) = {l/N){tT {z — v^)^^) are functional inverses of one another In the present 
context, we already noted that the ordinary complex variable z is replaced by the 
matrix Qq^ in ( [4.1| ), which leads us to consider matrix valued functions of Qq^ and 
composition thereof. Thus, following the nomenclature of the hermitean addition 
formalism, we name Ba the "Blue's" function associated with Q^'^l^. We now have 
to complete the arguments in the previous paragraph and show that composition of 
such matrix valued functions of matricial arguments are mathematically meaningful. 

In order that composition of such functions be a consistent operation (in the sense 
that it generalizes composition of complex- analytic functions of a complex variable), 
it is necessary that these matrix valued functions (with Qj^ in ( [4.4| ) being an important 
example), have the same matricial texture as Qq^ (Eq. (|4.1|) ). Recall that (Eq. ( [4.4|) ) 

/E 



— (rj — E)2 



V 



where rj and S = 2(77; r) are complex. Thus, for t] complex, the texture of Qj^ is more 
complicated than the texture of 

f Tj Z \ 



Go' 



\ z* r] J 



(Eq. ( |4.1| ), whose off-diagonal blocks are hermitean conjugates of each other. This 
mismatch of textures can be dealt with readily by generalizing the space of matrices 
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out of which Qq ^ is taken into 

^o"^ = I ^ I • (7.11) 

\ W 1] 

Here the off diagonal blocks are no longer related. To see that this generalization is 
the natural one to make, we follow the first few steps of analysis in Section 4 and find 
that with ( |7. 1 1| ) the self-energy matrix is 



Se = S(r/;z«;)5^ (7.12) 
(where the only change in ([4. 3]) is that now S depends on zw rather than on = z*z.) 



Thus, from the definition of self-energy ( [4.2|) we immediately see that ( [4.4|) is changed 
into 

1 f ^-v z \ ( g G 



g!^{v;z,w) 



zw — {rj — S)^ 



V w ^-V J \G' 



(7.13) 



This has essentialy the same texture as (|4.4|). With this extension, we have gained 



the fact that (|7.11|) , which generalizes here the ordinary complex variable z, has the 



same matricial texture as (|7.13| ). At the end of our calculations we will set w = z* 



We can also see that the new texture ( [7.11| ) is the right one at a more elementary 



level, by simply recognizing that the texture of ( [7.11| ) is preserved by raising this new 
g^^ to any integer (positive or negative) power. 

With this texture it is enough to consider and g^^ simply as 2 x 2 complex 
matrices with equal diagonal elements. Viewing ga {a = 1,2) as a 2 x 2 matrix 
valued function of a matricial argument ( [7.11| ) is an equivalent way of saying that 
we are dealing with the three functions ga{r]', zw), Ga{ri; z,w) and G'^{ri; z,w), and 
composition of such matrix valued functions amounts to compositions of such triads 
in the obvious way. 

In the following the symbols B, W, X etc. will stand for matrices of a texture 
identical to that of gQ^ and The matrix valued Blue's function Ba (W), associated 
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with the probabihty distribution Pa{4>l(t)a), (a = 1,2) is defined then as the matrix 

( Ba ba \ 

Ba=\ , (7.14) 

V K bJ 

that satisfies 

Ba[ga{X)]=X. (7.15) 

Similarly, we define the 2x2 matrix "Blue's" function B associated with ( |7. 5| ) by 
B [Q {X)] = X . From the definition of self-energy (|4.2| ), or equivalently from ( [7.7|) , 
we readily identify 

fi(>V) =>V-i + lE(>V) =>V-^ + l^r2J-trW) , (7.16) 

k=\ ^ 



where 1 is the 2x2 unit matrix. 

The identity of textures on both sides of ( [7.15| ) is assured already by the definitions 
of the quantities involved, and does not pose by itself any constraints on the various 
matrix elements that appear on both sides of ( |7.15| ). In other words, Ba is indeed the 
matrix valued inverse function of the matrix valued function Qa- Furthermore, due 
to this equality of textures we can readily apply Qa to both sides of (|7.15| ) obtaining 



Qa\Ba\Qa{Xy\\=Qa{X) , (7.17) 



which verifies that the matricial function composition in (|7.15| ) respects the obvious 



condition that taking the functional inverse twice in a row yields the original function. 
Thus, setting for example X = Qq^ in ( |7.15|) , that equation simply says that starting 



with the set of three functions QaiV'^ ^^)) G^a(^; z, w) and G'^{ri; z, w), we construct the 
set of three functions Ba{ri; zw),ba{ri; z,w) and b'^{ri; z,w), such that Qa{.Ba]b'J)a) = 
r], Ga{Ba] ba, &'„) = z and G'^{Ba\ ba, b'J = w. 

Note from ( |7l6| ) that i3(-W) = -B{W). This is so because B is the functional 
inverse of the propagator Q{X) which is manifestly an odd function of X due to the 
fact that in this paper all the probability distributions are even in 0. 
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Let us now go back to ( |7.8| ). The propagator on the right hand side of 

(|7.8|) is a function of ri,z and w, namely, of the matrix Qq^. We may general- 
ize and consider it a function of the variable W, rendering S^") a function of W, 

, , P -| 2k— 1 

namely, S('*)(W) = T,T=i'^2k [I tr G^^^W^ as in (|716|) . Similarly, by replacing 
the argument ^g^^ of on the right hand side of ( |7.9| ) by a variable A" (with X 
having the same texture as all the matrices above), we may consider the function 



2k-l 



Evaluating S(")(>V) at W = Wa = Ba [G {Go^)], and using ([7151 ) and we 
immediately find that 



2k 



E s("^(>Va) = E Er 

a k=l \ a / 

Thus, from ( [7.16| ) and from the matrix inverse of ( |7.15| ) we obtain 
E {Ba [G {X)] - G-' {X)\ = X- G-' {X) , 



(7.18) 



a=l,2 



namely. 



E Ba\G{X)\-G-\X) = X. 



(7.19) 



a=1.2 



Therefore, by definition, the Blue's function associated with (|7.5|) is simply 



Bt (W) = Bi (W) + B2 (W) - W"^ 



(7.20) 



in agreement with |^. 

To summarize, we carried out diagrammatic expansions (namely, expansions in 
traces of powers of G'^"'^^ and Gy^ of the self-energies S^^^ and S. We were allowed to 
do so because by construction G^^^ was taken to be non-hermitean, and thus in the 
analyticity domain of G^- Those diagrammatic expansions led us to introduce the 
matrix valued Blue's functions ( [7. 14] ) , in close analogy with the addition formalism 
for ordinary hermitean random matrices |jl^. As in the latter case, the Blue's func- 
tions follow a simple addition formula (|7.20|) under multiplication of the individual 
probability distributions (|7.5| ). 
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Now that we have Bt in ( [7.20| ), we can calculate the Green's function Q^{r]; z,w) 



associated with Ht by evaluating the functional inverse of Bt in (|7.2CI| ) at argument 



Gq^. After completing this step we may return to the original conventions and set 
w = z*. Also, for the sole purpose of determining the eigenvalue distribution of 
= 01 + 02 all we need is ^^(0 — iO; z, z*), and we may also set ?7 = — iO. 

7.1 A Simple Application: A Non-abelian Non- Hermit ean 
Central Limit Theorem 

Gauss proved that if we add K random numbers Xi, i = 1,2, ...K, with Xi taken from 
the probability distribution Pi{xi), then the normalized sum s = J2i follows the 
Gaussian distribution in the limit K tending to infinity. This result plays an important 
role in physics and mathematics and accounts for the ubiquitous appearance of the 
Gaussian distribution. 

A generalization of Gauss's central limit theorem for the case in which the variables 
Xi are replaced by x hermitean random matrices <fi, i = 1,2, ...K, taken from 
from the probability distributions 

P^{^^) = le-^*'-^^^^^). (7.21) 



certainly exists and was recently discussed in [jT4[-|3 -'■^ ^'^^^ version of the theorem 



the normalized sum = '^i indeed follows the Gaussian distribution in the 

large K limit. 

In this sub-section we wish to discuss Gauss's theorem for N x N non-hermitean 
random matrices (pi, i = 1,2, ...K, that are taken from from the probability distribu- 
tions 

P^{^U^) = ^e-^*'-^^(^I'^'). (7.22) 



^''Note that there exists another commonly considered class of random matrices, in which the 
element of the random matrices is each taken from a probability distribution (the same for each 
element). For this class, which is sometimes referred to as the Wigner distribution the theorem 
follows immediately from the usual abclian central limit theorem. Here we are speaking of the trace 



classes defined by (7.21) 
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As was already remarked in[0, it is not difficult to generalize one of the standard 
proofs of Gauss's theorem to matrices, and we shall do so at the end of this subsection. 
However, we find it quite interesting to address the question of whether (pT = Z^i 'Pi 
follows the Gaussian distribution by applying the non-hermitean addition formalism 
just presented. In doing so, we shall generalize the analysis made in the last section 
of [1^ to non-hermitean random matrices and show that Ginibre's disk distribution 



10| for the eigenvalues of 0t emerges naturally from the addition formalism. 

It is clear from the derivation of the addition formalism just given that when we 
add two random matrices, in general it is difficult to determine explicitly the resulting 

for the sum of the two matrices. What we hope for here is that the large K limit 
will bring considerable simplification. This is indeed the case. 

From the law of addition ( [7. 201) we learn that the function i3(VV) associated with 
the unknown Q{X) is given by 

B{W) = j:B^{W)-^—^. (7.23) 

4 = 1 

Using (|7.16|) we may rewrite (|7.23|) as 

oo / K \ / 1 \ 2fc-i 1 

^(w) = E Erg (otr>^) (7-24) 



.=1 v.=i / ^2 y w 

where is the 2k-th cumulant associated with the distribution Pj in ( |7.22| ). 

By definition (Eq. ( |7.14|) ), the Green's function Q{X) associated with Y.i4>i is 
determined by solving B{Q{X)) = X. We are however interested in the normalized 
sum 0r = Yl,i (pi- Define the corresponding propagator as 



QAX) ^ { {- 1 ) = VKQi^^X) , (7.25) 

where 
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Rescaling and substituting Q{y/KX) into ( |7.24| ), we find that Qt{'^) is determined 

by 

= X. (7.26) 



For finite K, ( [7. 261 ) is hopelessly complicated. However, as oo we see that 

it simplifies rather naturally to 

Qt\X) + ^ ^^Qt^X) = X (7.27) 



with 



r2 = ^Er?. (7.28) 



The self-energy associated with ( |7.271 ) is evidently = ^ tiQT{X). It is the sim- 

(\ 2k— 1 
I tr^j , where only 

r2 7^ 0. In other words, it is equal to the self-energy associated with the Gaussian 
distribution 

Pg(0V) = ^exp (^-^ tr0t0) . (7.29) 



Re-substituting 



V 



and using the general hermitization formalism in Section 2 (and in particular Eqs. 
(|373|) and ( p.l2[ )) we conclude |@, |^ that the eigenvalues of 0t = (1/V^) Y.f 4>i are 
distributed uniformly in a disk of radius \/T2. 

Thus, not only do we obtain Ginibre's distribution, we also learn that r2 is given 
by (HH). 



While we have proved that the density of eigenvalues of (^t is identical to the 
density associated with the Gaussian distribution (|7.29|) , we cannot yet conclude 



that the probability distribution of (f)T is Gaussian. The reason is that in principle, 
a given eigenvalue distribution may correspond to more than a single probability 
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distribution^. To see that is indeed taken from a Gaussian distribution, we have 
to calculate that distribution explicitly. 

As mentioned above, it is not difficult to extend one of the usual proofs of the 
central limit theorem to the case of matrices. The distribution for the normalized 
sum matrix = X^i 0i is given by (here we omit the subscript T) 

K 

I dd.PAd\d.) \5\d--, 

'K 



= / (n/ d<P^Pi<PU^)j e^^^ '"^ "t' . (7.30) 

The integral over the (pi can be done in the large K limit: 

= 1 - ^ / d<l>,Pi4<l>.) tr m tr(tt0t) + o (-i^) 



J- 2 

4KN 



tr (th) + O (-^) (7.31) 



where by definition r2 = / d(f)iP{(f)l(j)i)jj tr^f^j. Reexponentiating and integrating 
over t we obtain the desired result that P{(f)) is proportional to e ^ tr</>t<}i^ with r2 
given by (|7:28|) . 
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^'''This is indeed the case for hermitean random matrices. For example, for matrices in the Wigner 
class, even when the distribution is not Gaussian, the corresponding density of eigenvalues still 
satisfies the semi-circle law, as is well known. 
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Appendix : A Note Concerning the "Single Ring Theorem" 

The hermitean matrix {/)^(/) can always be diagonahzed = h?U by a unitary 
matrix [/, with = diag(A^, A|, • • • , A^), where the Aj are all real. This implies 
that = V^^At/, with V a unitary matrix as well. Thus, the complex eigenvalues 
of (f) are given by the roots of det(z — NW) = 0, with W = UV'^ . Evidently, as 
W ranges over U{N), the eigenvalues of AW could be smeared (in the sense that 
they would not span narrow annuli around the circles |2;| = |Aj|.) In Section 6, we 
proved (in the large N limit) that for taken from the class of ensembles of the form 
P{(f)^(f)) = (1/Z)exp[— A^tr V{(f)''(j))] (Eq. ( |1.1D ), the eigenvalues can only span either 



a disk or an annulus.We have verified that this is the case, explicitly for N = 2, and 
numerically for large A^ (of the order of 100) and for some range of A^. 
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